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Abstract 

We explore the cosmological evolution of equation of state (EoS) for dark energy in g-essence 
f) i' models, the action of which is described by a function of both the canonical kinetic term of both 



the scalar and fermionic fields. We examine g-essence models with realizing the crossing of the 
phantom divide line of ii^de = — 1 as well as the models in which the universe always stays in the 
non-phantom (quintessence) phase (tuoE > ~1)- In particular, we find an explicit g-essence model 
with the crossing from the non-phantom phase to the phantom one {wde < — 1). This transition 
behavior is consistent with the recent observational data analyses. 
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I. INTRODUCTION 



, 2|, cosmic 



According to recent cosmological observations such a. Type la Supernovae Q 

microwave background (CMB) radiation |3|-[7|, large scale structure (LSS) 8|, |9|, baryon 
acoustic oscillations (BAO) 10|, and weak lensing ll|, the expansion of the current universe 
is accelerating. 

There are two categories to explain the current cosmic acceleration. One is to suppose 
the existence of some unknown matters called dark energy. The other is to modify gravity. 
There are a number of candidates for dark ener gy, e.g., the cosmological constant, scalar 



fields like quintessence |l2-|l6| and k-essence |17H22|. and fluid as a Chaplygin gas (23|, l24| 
(for reviews, see, e.g., J25l-l32| . and for related works, see, for example, SSHSTJ). There also 
exist several ways of modification of gravity, e.g., F{R) gravity (for recent reviews, see, 



e.g. 



38 



45|). 



One of the most important quantities in the studies in terms of the property of dark 
energy is the equation of state (EoS). It seem to be suggested by recent various cosmological 
and astronomical observational data 46N51| that the crossing of the phantom divide line of 
wde = — 1 occurred in the near past, where wde is the EoS for dark energy. There are a lot 
of theoretical attempts and proposals to account for such a phantom crossing phenomenon 
in the framework of both dark energy in general relativity and modified gravity scenarios in 



25 



31 



38 



m- 



the literature (for references, see, e.g., the review articles shown above 

In addition, on the theoretical studies, as phenomenological generalizations of k-essence 
models, the action of which is written by a function of the canonical kinetic term of a scalar 
field, f-essence and g-essence models have recently been proposed in Ref. 52|. The action 
of f-essence is described by a function of the canonical kinetic term of a fermionic field. 
On the other hand, the action of g-essence is represented by a function of the canonical 
kinetic term of a scalar filed as well as that of a fermionic field, which corresponds to a 
generalized theory constructed by combining k-essence with f-essence. In this paper, we 
concentrate on the evolutions of cosmological quantities, especially the EoS for dark energy 
in g-essence models. We analyze the evolutions of cosmological quantities in g-essence models 
with realizing the crossing of the phantom divide as well as without it. In particular, we 
examine g-essence models in which the universe always evolves within the non-phantom 
(quintessence) phase (wde > ~1)- In addition, we find an explicit g-essence model with 



realizing the crossing of the phantom divide hne of wde = — 1 from the non-phantom phase 
to the phantom one (u'de < — !)• It is interesting to emphasize that this transition behavior 



46 



It has 



is compatible with the data analyses of the recent cosmological observations 

also been investigated that in viable f{R) gravity models, the phantom crossing with the 

opposite manner (i.e., from the phantom phase to the non-phantom phantom one) can occur 



(see, e.g. 



53| and the references therein), whereas in f(T) theories, where T is the torsion 



scalar, the phantom crossing with the sarne transition direction (i.e., from the non-phantom 
phase to the phantom one) can happen 

e.g., yjsgl). 



siise 



56| (for related aspects in /(T) theories, see. 



The most significant physical motivation of this work is to illustrate the cosmological 
evolution of the EoS for dark energy wde in g-essence models. This is inspired by the data 
analysis of recent cosmological and astronomical observations 46l-l5l|. which implies that the 
crossing of the phantom divide happened in the near past. On the other hand, it is known 



that in the so-called "reconstruction" 



60 



69|. we starts with a well motivated theoretical 
model such as simplicity and/or analogy with some fundamental theory, and by using its 
form derives an expansion law of the universe to be compared with observational data, or 
inversely, starts with the observational data themselves and try to fit them with a theoretical 
model. Our demonstrations in this work are complementary to such reconstruction processes. 
The paper is organized as follows. In Sec. II, we briefiy explain g-essence, k-essence and 
f-essence models. In Sec. Ill, we investigate a particular g-essence model and the evolutions 
of cosmological quantities. In Sec. IV, we explore g-essence models both with and without 
the crossing of the phantom divide. Finally, conclusions are presented in Sec. V. 



II. ESSENCE MODELS 



In this section, we explain g-essence, k-essence and f-essence models. 



A. g-essence 



The action of g-essence is given by 52 1 



S= d^x^^[R + 2K{X, r, (j), ^, ^)], 



(2.1) 



where K is some function of its arguments, is a scalar function, if) = {ilji,ip2,'ip3,'ip4:)'^ is a 
fermionic function and ip = ip'^'-f^ is its adjoint function.® Here, the canonical kinetic term 
for the scalar field X and fermionic field Y are given by 



X = O.S^^^'^V^^V,^ , Y = 0.5t[ijT^^D^^ - iD^ij)Ty] 



(2.2) 



where V^ and D^ are the covariant derivatives. We note that here the fermionic fields are 
treated as classically commuting fields. 

We consider the fiat Friedmann-Lemaitre- Robertson- Walker (FLRW) space-time with the 
metric, 

ds^ = dt^ - a^{dx^ + dy^ + dz^) , (2.3) 



where a{t) is the scale factor, and the vierbein is chosen to be (see, e.g., |70| ) 



«) = diag(l, 1/a, 1/a, 1/a) , (e°) = diag(l, a, a, a) . 



(2.4) 



In the case of the fiat FLRW metric (2.3), from the action (2.1) the basic equations 



read 



52| 



3H^ -p = , 

2H + 3H^+p = , 

Kx4> + {Kx + ^HKx)^ -K^ = 0, 

Ky^ + 0.5{3HKy + KY)iJ - i-f'^K^ = 0, 

Ky'$ + 0.5{3HKy + kY)ip + iK^-f^ = 0, 

p + 3H{p + p) = 0, 

where the kinetic terms, the energy density and pressure take the forms 



(2.5) 
(2.6) 
(2.7) 
(2.8) 
(2.9) 
(2.10) 



X = 0.5</)2, F = 0.5z(^7°^ - ^7V) , p = 2KxX + KyY - K , p = K. (2.11) 

Here, a dot denotes a time derivative of d/dt and H = a/a is the Hubble parameter. 
Moreover, the subscription of K, e.g., Kx = dK/dX, denotes the derivative of K with 
respect to X. Several properties of g-essence have been studied in Refs. 7ll-l74l| . We remark 
that the model (2.1) can also describe k-essence and f-essence, as is shown in the following 
next subsections. 



B. k-essence 



We examine the case that the Lagrangian K has the form 



K = Ki = Ki{X,> 



(2.12) 



Then, the action (2.1) takes the form 52 1 



S= f d^Xy/^[R + 2Ki{X,, 



(2.13) 



which corresponds to k-essence 17H22| (for a solvable case, see 75|). For the flat FLRW 
metric (2.3), the basic equations of k-essence becomes |l7H22| 



3H'-p 



2H + 3H'+p 
p + 3H{p + p) 



0, 
0, 
0, 
0. 



(2.14) 
(2.15) 
(2.16) 
(2.17) 



C. f-essence 



Now, we consider the case that the Lagrangian K takes the form 



K = K2 = K2{Y,ij,2P) 



In this case, the action (2.1) is written as 



S = / d^xy^[R + 2K2{Y, iP, ij)] 



(2.18) 



(2.19) 



This is the so-called f-essence 52|]. In the case of the flat FLRW metric (2.3), the basic 
equations in f-essence have the forms [52 1 



3H'-p 



2H + 3H'+p 
Kvi} + 0.5{3HKy + Ky)iIj - il^K^ 
Kyi + 0.5{3HKy + Ky)^ + iK^a^ 

p + 3H{p + p) 











(2.20) 
(2.21) 
(2.22) 
(2.23) 
(2.24) 



where the energy density and pressure are expressed as 

P = KyY-K, p = K. (2.25) 

It is important to note that in addition to the description of a spinor field of f-essence as 
a classical c-number quantity, it can also be treated as a Grassmann-valued quantity (see, 



e.g. 



76| ) , or as an operator ( "q-number" ) , as in Ref . 77 1 . 



III. PARTICULAR G-ESSENCE MODEL 



In general, in g-essence models 52|] equations of motion are very complicated and therefore 
it is difficult to obtain those exact analytical solutions. In contrast to such more general and 
complex cases, in order to execute the analyses we here explore the following more simple 
particular model 

K = eX + (rY-Vi{^)-V2{u), (3.1) 

where u = ipip, and e and a are constants. In this model, the equation system (2.9)-(2.14) 
is described as 

3H^- p = 0, (3.2) 

2H + 3H'^+p = 0, (3.3) 

e4> + 3eH^ + Vi^ = 0, (3.4) 

a^ + -aiJV^ + z^aW = 0, (3.5) 

aip + -aHip - iV^ip-r^ = , (3.6) 

p + 3H{p + p) = 0, (3.7) 

where 

p = O.5e02 + Vi + V^2 , (3.8) 

p = O.5e0^ - 14 - 14 + V^u . (3.9) 

Here, a prime denotes the derivative of a quantity with respect of its argument such as 
V = dV/du. We note that the contribution of the fermion field is included in the potential 

V2{U). 



IV. COSMOLOGICAL SOLUTIONS 

In this section, we analyze cosmological solutions in several examples with Vi = in 
the particular g-essence model shown in Sec. III. In particular, we examine the cosmological 
evolutions of the energy density p and pressure p, the deceleration parameter q = —a/ {aH"^), 



the jerk parameter j = 'a / (aH^) [78|, and the EoS w = p/p as functions of t. We study both 
four examples with and those without the crossing of the phantom divide. We note that 
the subscription i = 1, 2, . . . , 8 oi pi, pi, qi, ji and Wi denotes the number of four examples 
in Sees. IV A 1-4 (z = l,2,3,4)and those in IV B 1-4 (z = 5,6,7,8). In the A Cold Dark 
Matter (CDM) model, q = —l,j = 1 and w = —1, and hence these quantities denote the 
deviation of a model from the ACDM model. In what follows, we consider the dark energy 
dominated stage in which dark energy (which is, e.g., g-essence in this paper) is dominant 
over the universe, and therefore it can be interpreted that w ~ wde, where wbe is the EoS 
for dark energy. 

It is known that in the FLRW background (2.3), the effective EoS for the universe 
is given by |38|, l39| Wcs = Pcs/PcS = -1 - 2H/ {3H'^) with pcfr = SH^/k"^ and p^s = 



2H + 3H^ j /k^ being the effective energy density and pressure of the universe. Here, 
Pcff and Pcs correspond to the total energy density and pressure of the universe, respectively. 
If the energy density of dark energy is dominant over that of matter completely, wbe ~ Wcs- 
In the non-phantom (quintessence) phase, H < and hence Wcs > —1, whereas in the 
phantom phase, H > and therefore Wcs < —1. Moreover, for H = 0, Wcs = —1, which 
corresponds to the cosmological constant. From the above considerations, we examine the 
case that w ~ wde ~ Wcs- 

In case of Vi = 0, the system (I3.2l) - (l3.7p becomes 

3H^-p = 0, (4.1) 

2H + 3H^+p = 0, (4.2) 

4) + 3H(j) = 0, (4.3) 

atp + -aHi! + iV^'y^ilj = 0, (4.4) 

aij + -aHtP - iV'^^-f = , (4.5) 

p + 3H{p + p) = 0, (4.6) 



where 

p = O.5e02 + V2 , (4.7) 

p = O.5e02 -V2 + V^u . (4.8) 

As a result, we find 

„=4, j,^tl, ^,._^^.3ff^ (4.9) 

and 

V^/ = ^e^^, ^. = ^e-^, (4.10) 

where c and 0o are constants. Here, / = 0, 1, /c = 2, 3, cj (j = 0, . . . , 3) obey the following 
condition c = |cop + |cip — |c2p — |c3p, and 

D = — (2 I a^dH + (pl I a-'^dtj . (4.11) 

A. Solutions without the crossing of the phantom divide 

In this subsection, we explore derive the solutions in four examples without the crossing 
of the phantom divide. 

1. Example 1 

We start from the following expression for the scale factor 

2 dt 

■I 



a = e 



where 



3 ' t - 0.5(t2 sin[f ] + t cos[i] + Si[^]) .^ ^2) 



Si[x]= / ^-^^dz. (4.13) 

Jo ^ 



In this case, we obtain 



H = (A 14) 

3[t-0.5(t2sin[i] + tcos[i] + S2[i])] ' ^' ' 

dt 

^ = ,i"^t-0.5{t^sm[^+tcos[^ + St{^]) ^ (4.15) 

dt 

= 0oe''^^-O-^(^'^i^[l]+^^°^[l] + ^^[l]), (4.16) 

dt 
^, = Qe^''"^^-0-5(^^sin[l]+tcos[l] + 5^[i]) ^^_^^^ 

dt 

^, = ce"^"""^ t-0.5(t2sin[l]+tcos[i] + 5^[l]) ^ ^^^^g^ 

where 

^^^ /"r 4(l-tsin[|]) '^ t - 0.5(t2 sin[i] + t cos[|] + Si[\]) 

ca J ^3[t-0.5{t^sm[\] + tcos[\] + Si[]])]^ 

dt 



(4.19) 

The corresponding potential has the form 

dt 

V, = -O.5e0ge''^^ - 0-5(^2 sin[f] + tcos[i] + 5z[f]) 

4 
^ 3[t - 0.5(t2 sin[i] + t cos[i] + 5z[i])]2 ' ^^'^^^ 

The energy density and pressure are described by 

4 
Pi 



3[t - 0.5(t2 sin[i] + t cos[i] + St[^])]^ ' 
At sin i 
^' " ~3[t - 0.5(t2sin[i] + tcos[l] + Si[^W ' 

For this example, the EoS and the deceleration parameter have the form 

W7i = — = -tsin[-], (4.22) 

Pi t 

and 

q, = ^ = 0.5 - 1.5t sin[i] . (4.23) 

all t 




FIG. 1: pi (solid line) and pi (dotted line) as a functions of t. 
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FIG. 2: wi as a function of t. 

In Fig. [H we show the cosmological evolutions of the energy density pi and pressure pi as 
the functions of t. Furthermore, in Fig. [2] we depict the cosmological evolutions of the EoS 
Wi as a function of t. From Fig. |2] we see that in this model the crossing of the phantom 
divide cannot be realized. In addition, the jerk parameter is defined as 

i2 



Ji = 



1-1 

8t 



a a 



■ L 



rl. 



{2 - 3r) sin'l-] + Si[-]{t sm[-] - cos[-]) + tcos[-](2 - cos[-]) 

V L L (j L L 



(4.24) 



In Fig. [3l we plot the deceleration parameter qi and jerk parameter ji as functions of t. 



2. Example 2 



We suppose that the scale factor is described by the following expression 



2 



dt 



a = e 



3' t-eO-5*"i(t + l)0-5(*'-i)i-o- 



5*2 



(4.25) 
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FIG. 3: qi (solid line) and j'l (dotted line) as functions of t. 



Then, we find 



H 



u = c e 



3[t-eO-5*-i(i+l)0-5(t^-i)t-o.5<2]' 
dt 

~^^t- e0-5*-l(t + I)0.5(t2-l)^-0.5t2 



-2/ 



di 



5oe 



tpk 



iD-J 



cie 



Cfce 



^_g0.5t-l(^+l)0.5(t2-l)^-0.5i2 
i_e0.5t-l(^+l)0.5(t2-l)^-0.5t2 



iD-/ 



dt 



t-eO-5*-i(t + l)0-5(*'-i)t-o- 



5t2 



(4.26) 

(4.27) 
(4.28) 
(4.29) 
(4.30) 



where 



D 



^ 4(_t _ 1 + e0-5t-i(t + i)0.5{t2+i)^-o.5t2 In -^ + ^ ' 



C 0" 



2/ 



t 



X e 



3(tv/^TT - e0-5*-i(t + 1)0.5*2^-0.5*2) 

dt dt 

t - e0-5*-l(t + I)0.5{t2-l)^-0.5t2 ^ ^2 f Q~^^t- e0-5*-l(t + l)0-5(t^-l)t"0.5t2 ^^^ ^^^^ 



The form of the corresponding potential is given by 

dt 

-4/: 



r. r l2 -^ + _ p0.5i-l('+ I 1^0.5(t2-l) + -0.5t2 



3[i_e0.5i-l(^+l)0.5(i2_l)^»0.5t2p 



r- (4.32) 



The energy density and pressure are expressed as 

4 



P2 



3(i_e0.5t-l(i+l)0.5(t2-l)^_0.5i2)2' 



P2 



4(1 + i)* 



3e(t - e0-5*-i(^ + l)0-5(*'"i)t-o-5*')2 ' 

(4.33) 



11 




FIG. 4: p2 (solid line) and p2 (dotted line) as functions of t. 




FIG. 5: W2 as a, function of t. 



In case of this example, the EoS and the deceleration parameter become 



W2 



'Hi 



(4.34) 



and 



g2 = 0.5-1.5e~^(l 



(4.35) 



In Fig. m we illustrate the cosmological evolutions of the energy density p2 and pressure 
P2 as functions of t. In Fig. |5l we show the cosmological evolutions of the EoS W2 as a 
function of t. It follows from Fig. [5] that the crossing of the phantom divide cannot occur 
in this example. In Fig. |6l we plot the deceleration parameter q2 and jerk parameter J2^ as 
functions of t. 



^ In Fig.ini we have plotted the behavior of J2 numerically because the analytic expression is too complicated. 
This is the same as for J3 in Fig. [TU]in Sec. IV A 3. 



12 




FIG. 6: 52 (solid line) and J2 (dotted line) as functions of t. 



3. Example 3 



We provide that the scale factor is expressed as 



dt 



a = e 



3 ■' ^ _ j\^0.5A^QA{0.5t-l)U _|_ y^-\0.5{t2-A2)^-0.5t2 



(4.36) 



where A is a constant. From this expression, we have 
H ' 



dt 



-2/ 



dt 



5oe 



f _ j\0.5A^ QA{0.5t-l) U _|_ j^\0.5{t^-A^)-f.-0.5t^ 

dt 



^^ ^ ^^^^ ^t_^0.5A2gA(0.5i-l)(i + ^)0.5(i2-A2)^-0.5i2 

dt 

^^ ^ ^^g"'^"^ t - A0-5^'e^(0-5*-l)(t + yl)0.5(t2-A2)^-0.5t2 



(4.37) 

(4.38) 
(4.39) 
(4.40) 
(4.41) 



where 



C (J 



4(_1 -|_ ^0.5^2 gA(0.5i-l)(^^ _^ ^)0.5(t2_A2)^-0.5i2+l Jn .^ "^ "^ ^ 



t 



3('^ _ ^0.5A2gA(0.5i-l)C^ _|_ ^~\0.5(t2-A2)^-0.5t2N2 



+00 e 



-2/ 



rft 



t - A0-5^'e^(0-5*-i)(t + A)0-5(*^-^^)t-o-5*^ 1^^ _ ^4 42) 
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FIG. 7: pz (solid line) and pa (dotted line) as functions of t for A = 10. 



The corresponding potential becomes 



dt 



V2 = -0.5e</.2e '^ ^ - AO-5^'e^(o-5*-i)(t + A)0-5(*'-^')t-o-5*' ^ 

4 

3[t — yl0.5A2gA(0.5t-l)C^ _|_ ^^0.5(t2-A2)^-0.5i2l2 ' 

The energy density and pressure are given by 

4 



P3 



P3 



3[^ _ ^0.5A2gA(0.5t-l)C^ _|_ ^^0.5(t2-A2)^-0.5t2l5 

4(1 + 4)* 



3eAr^ _ ^0.5A2gA(0.5t-l)C^ _|_ ^^0.5(t2-A2)^-0.5t2l2 ■ 

In this example, the EoS and the deceleration parameter are written as 



W3 



;i + y)' 



and 



gg = 0.5 - 1.5e-^(l + 



A. 



(4.43) 



(4.44) 
(4.45) 



(4.46) 



(4.47) 



respectively. 

In Fig. [TJ we plot the cosmological evolutions of the energy density ps and pressure ps as 
functions of t for A = 10. In Fig. |8l we display the cosmological evolution of the EoS W3 as a 
function of t for A = 10. We find from Fig. [8] that in this model the crossing of the phantom 
divide cannot happen. In Fig. [9l we show the cosmological evolution of the deceleration 
parameter q^ as a function of t. In Fig. [TOl we illustrate the cosmological evolution of the 
jerk parameter J3^ as a function of t. 



2 See the footnote in Sec. IV A 2. 
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FIG. 8: W3 as a function of t for A = 10. 




FIG. 9: qs (solid line) as a function of t. 



4- Example 4 



We consider that the scale factor is written by 



a = {ao + b cosh [kt] + d sinh [kt] ) 3 . 



(4.48) 



where Oq, b, d and k are constants. In this case, we acquire 

k{b siVihlkt] + dcos\i[kt]) 



H 



u 



3(ao + 6 cosh [H] + (isinh[/ct]) 
c 

ao + 6 cosh [H] + (isinh[fct] 

2(1)1 ^.,„,^„ ( (p- ao) tanh[0.5A;t] + d 



tpk 



Cl 



: arctan 



v/&2 - d2 



JD 



(oo + 6 cosh [H] + (isinh[H])2 

Cfc 

(flo + 6cosh[A;t] + (isinh[A;t])^ 



iD 



(4.49) 
(4.50) 

(4.51) 

(4.52) 
(4.53) 



15 



50 




40 




so- 


/ 


la 


■ 


20 
10 


/ 







12 3 4 



FIG. 10: js (dotted line) as a function of t. 



where 



D 



2k 
3c a 



nb^-d'--l + ^. 






^/b^ - d? 



arctan 



{h — ao) tanh[0.5A;t] + d 



Vb' 



'2 - rf2 _ al 



/tanh[0.5H+l] 
+ "°^''Unh[0.5H-l] 



The corresponding potential is given by 

V2 = 



O.5e0g 



fc2 / 6sinh[H] +(icosh[H] V 
ao + b cosh [kt] + d sinh [kt]y 3 \ Oq + 6 cosh [kt] + d sinh [kt] J 



The energy density and pressure have the forms 



P4 
P4 



k'^ f 6sinh[/i;t] + (icosh[A;t] 



3 V'^o + bcosh[kt] + (isinh[fct] 
P (6 sinh [H] +dcosh[kt])'^ - 2[ao(6cosh[A;t] + dsmh[kt\) + b"^ - d"^ 



3 (ao + 6 cosh [H] + dsmh.[kt]y 

For this example, the EoS and the deceleration parameter become 

2[ao(6cosh[H] +dsmh[kt]) + b"^ - d'^] 



and 



004 



Qi 



4 



-1 



(6 sinh [H] +cicosh[H])2 
3[ao(6cosh[A;t] + dsinh[H]) + b^ - d"^] 



(4.54) 



(4.55) 



(4.56) 
(4.57) 



(4.58) 



(4.59) 



(6sinh[H] +rfcosh[/ct])2 
respectively. 

In Fig. [HI we depict the cosmological evolutions of the energy density p^ and pressure p^ 
as functions of t. In Fig. [121 we show the cosmological evolution of the EoS W4 as a function 
of t. From Fig. [121 we understand that also in this example the crossing of the phantom 
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FIG. 11: p4 (solid line) and p4 (dotted line) as functions of t. 




(4.60) 



FIG. 12: W4 as a function of t for uq = —1/2, b = 1/2, d = 1 and k = \/3. 

divide cannot be realized. Moreover, the jerk parameter is described by 

_9al- 10^2 + 9rf2 + (^2 + t/2) cosh2[H] + bdsinh[2kt] 
^^~ (6sinh[H]+rfcosh[H])2 ' 

In Fig. [131 we illustrate the cosmological evolutions of the deceleration parameter 94 and 
jerk parameter J4 as functions of t. 

In summary, from Figs. [21 [5l [8] and [12] we see that in the above four models, the universe 
is always in the non-phantom (quintessence) phase (in which w > —1) and therefore the 
crossing of the phantom divide cannot be realized. 



B. Solutions with the crossing of the phantom divide 

Next, in this subsection we investigate the solutions in four examples with realizing the 
crossing of the phantom divide. 
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FIG. 13: q4 (solid line) and J4 (dotted line) as functions of t. 



1. Example 1 



We consider that the scale factor is described by 

dt 



a = e 



3^ t- Jttan[\]dt 



(4.61) 



In this example, we acquire 



H 



3(t-/ttan[i]t/t) 



-2/ 



dt 



u = c e 



-2/ 



5oe 



iD-f 






cie 



Cke 



t- Jttan[\]dt 

dt 
t- Jttsin[\]dt 

dt 
t - f ttan[]]dt 



^D^I 



dt 



t - j tta.\\[\]dt 



(4.62) 

(4.63) 
(4.64) 
(4.65) 
(4.66) 



where 



D 



c a 



4(ttan[i] -1) 
3(t-/ttan[i]dt)^ 



2/ 



dt 



2/' 



dt 



^ I 

t — / ttan[j](it I i2 t — j ttaxi[j\dt 



dt. (4.67) 



The form of the corresponding potential reads 

dt 

4/- 



V2 = -Q. feeble ^-/^tan[i]rft 



tJ--^ + 



?>{t- jti&n[\]dty 



(4.68) 
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FIG. 14: It's as a function of t. 



The energy density and pressure are given by 

4 

P5 



4ttan[il 



3(t-/ttan[i]rft)2' 



P5 



3(t-/ttan[i]dt)^ 



For this example, the EoS and the deceleration parameter are written by 



and 



W5 = — ttan[-] 



gs = 0.5- 1.5ttan[-] 



(4.69) 



(4.70) 



(4.71) 



respectively. 

In Fig. [TH we display the EoS w^^ as a function of t. From Fig. [TU we see that in this 
model the crossing of the phantom divide can be realized. In addition, the jerk parameter 
is described by 

J5 = 5t + 9ttan[-] (t- (2t^ - l)tan[-] j +9 (ttan[-] -tan^[-] - 1 J / ttan[-]rft . 

(4.72) 



2. Example 2 



We take the scale factor as 



a = e 



2 dt 

3^t-0.5(t2sin[i] + tcos[i]) 



(4.73) 
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In this case, we find 



H = (A 7A) 

3[t-0.5(t2sin[i] + tcos[^])] ' ^ ' ' 

dt 
u = c e"'^t-0.5(t2sin[l]+tcos[i])^ ^^^^^ 

dt 

dt 

^, = cj''~^ t-0.5(t2sin[l]+tcos[i]) (477) 

dt 
^, = ce"^"""^ t-0.5(t2sin[i]+tcos[i]) ^ ^^^g^ 



where 



^^^ /• 4(l-[t + 0.5t-^]sin[|]) '^t-0.5(t2sin[l]+tcos[i]) , 

cai ^3(t-0.5(t2sin[i] + tcos[i]))2 "" 

-2/ 



t-0.5(t2sin[i] + tcos[i]) 



+0^ / e '- "-^^^ ""H t J "^ ' ^^"^t J^ jdt . (4.79) 



The corresponding potential is given by 

dt 

-4/- 



V, = -O.5e02e~ ^ " 0-5(^' sin[i] + t cos[i]) ^ ^ ^^gQ^ 

^° 3[t-0.5(t2sin[i] + tcos[i])]2 ^ ' 

The forms of the energy density and pressure become 

^' 3[t-0.5(t2sin[i]+tcos[i])]2' ^' 3[t - 0.5(^2 sin[i] + tcos[i])]2 ' ^' ^ 

For this example, the EoS and the deceleration parameter are given by 

We = -{t + O.Sr^) sin[-] (4.82) 



and 



qe = 0.5 - 1.5(t + 0.5t"^) sin[-] , (4.83) 



respectively. 

In Fig. [15], we display the cosmological evolution of the energy density pe and pressure pe 
as functions of t. In Fig. [11], we show the cosmological evolution of the EoS Wq as a function 
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FIG. 15: pq (solid line) and pg (dotted line) as functions of t. 




FIG. 16: vjQ as a function of t. 

of t. It follows from Fig. [16] that in this case the crossing of the phantom divide can occur. 
Furthermore, the jerk parameter is expressed as 



J6 = Y^[18 + 97t' + 9t sin[^] + 18^(3^^ sin[^] - 21" - 3) sin[^] 

-9(cos[-](7t2 + 1) + At^ + 2) cos[-] 



(4.84) 



In Fig. [T71 we demonstrate the cosmological evolutions of the deceleration parameter gg and 
jerk parameter jg as functions of t. 

3. Example 3 

We express the scale factor as 



a = flo + &sin[/i;t] + dcos[kt] 



(4.85) 
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FIG. 17: qq (solid line) and Jq (dotted line) as functions of t 

where ao, fc, d and k are constants. From this description, we obtain 

k{bcos[kt] ~dsm[kt]) 



H 



u 



^i 



i^k 



ao + 6sin[A;t] + dcos[kt] 
c 

(ao + 6sin[A;t] + dcos[kt])^ ' 

00 

(ao + 6sin[A;t] + dcos[kt])^ ' 

Q 

(ao + 6sin[/i;t] + dcos[kt]y-^ 

Ck 

(ao + 6sin[A;t] + dcos[kt]y-^ 



JD 



-iD 



where 



D = — / (2A;^(ao + 6sin[H] + rfcos[A;t])(ao(6sin[H] +dcos[kt]) + 6^ + rf^) + 



The corresponding potential is given by 



Vo 



O.5e02 



+ 3fc^ 



6cos[A;t] — (isin[/i;t] 



(ao + 6sin[H] + dco?,[kt\Y \ao + 6sin[A;t] + dcos[kt\ 

The energy density and the pressure have the forms 



P7 = 3A;2 
p-j = -k^ 



6cos[A;t] — (isin[/i;t] 
aQ + b sm[kt] + d cos[kt] 
3(6 cos [H] - dsm[kt]f - 2[ao(fesin[H] +dcos[kt]) + 6^ + d^] 



(ao + 6sin[A;t] + dcos[kt]y 
In this example, the EoS and the deceleration parameter are described as 



_ 2 [ao (6 sin [kt] + d cos [kt] ) +b'^ + d^ 

^^ ~ ~ ^ 3(6cos[H]-rfsin[A;t])2 



(4.86) 

(4.87) 

(4.88) 
(4.89) 
(4.90) 



(ao + fesin[H] + dcos[kt])-^)dt . (4.91) 



(4.92) 



(4.93) 
(4.94) 

(4.95) 
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FIG. 18: p7 (solid line) and p^ (dotted line) as a functions of t for uq = 2, b = 1, d = 1 and k = 2. 
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FIG. 19: wj as a function of t for oq = 2, 6 =, d = 1 and k = 2. 



and 



g? 



-1 



[ao(6sin[H] -F rfcos[/ct]) -h 6^ + rf^ 



(4.96) 



{b cos [/ct] — d sin [/ct] ) ^ 
respectively. 

In Fig. [iHl we show the cosmological evolutions of the energy density pj and pressure p^ 
as the functions of t. In Fig. [191 we illustrate the cosmological evolution of the EoS wj as a 
function of t. We find from Fig. [19] that in this model the crossing of the phantom divide 
can happen. Furthermore, the jerk parameter is written by 

(oq + 6sin[/i;t] -|- dcoslkt])"^ 



J7 



(4.97) 



(62 - d^) cos^[kt] - bdsm[2kt] + d^ ■ 

In Fig. [201 we display the cosmological evolutions of the deceleration parameter qj and jerk 
parameter J7 as functions of t. 
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FIG. 20: q2 (solid line) and J2 (dotted line) as functions of t. 



4- Example 4 



We start from the following expression for the scale factor 

2 dt 



(4.98) 



where / = Oq + bcosh.[kt] + (isinh[/i;t], and /, m, n, Qq, b, d and k are constants. Then, we 
obtain 



H 



3[t + f{lf-^ + mf~^ + n)dt] ' 
dt 



2/: 



rft 



dt 
^, = cie t + j{lt^ + mt^ + n)dt 



iD-J 



dt 



^, = c,e 't + J{ir- + mr^ + n)dt 



(4.99) 

(4.100) 
(4.101) 
(4.102) 
(4.103) 



where 



D = — 

c a 



dt 
4(l + //-2 + m/-^+n) ^■^t + J{lf-'' + mf-^ + n)dt, 

3it + J{lf-^ + mf-^ + n)dt)^ 

dt 



(4.104) 
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The corresponding potential becomes 

dt 



v.. . -,^.^'^ * + m-' + "/-'+ n)M ^ 3,^^j^,^,,;,„^_.^„^^^,, (4.105) 

The energy density and pressure are given by 

4 

^' " 3[t + /(//"2 + ^j-i+^)rf^]2 ' (4.106) 

In this example, the EoS and the deceleration parameter have the following forms 

u;8 = //-2 + ^j-i + ^ (4.108) 

and 

g8 = 0.5 + 1.5(//-2 + m/"i + n), (4.109) 

or 

w^ = /(flo + bcosh[kt] + rfsinh[/ct])~^ + m(ao + 6cosh[H] + rfsinh[H])~^ + n (4.110) 

and 

gg = 0.5 + 1.5(/(ao + &cosh[H] + t/sinh[H])"^ + m(ao + &cosh[H]+c?sinh[H])"^+n). (4.111) 

In Fig. [21], we illustrate the EoS wg as a function of t. From Fig. [21], we understand 
that in this case the crossing of the phantom divide can be realized from the non-phantom 
(quintessence) phase (in which ty > — 1) to the phantom one (in which w < —1). Moreover, 
the jerk parameter is described by 

is = JTzWfl + fm + /2) + 36//(m + n) + 18m f{m + 2fn) + 2/^(2 + 9n^) + 

+9ff(2l + fm)[t+ f(lf-^ + mr' + n)dt]]. (4.112) 



As a consequence, from Figs. [2] [S] [H] and [12] we understand that in the above four models, 
the universe is always in the non-phantom (quintessence) phase and therefore the crossing 
of the phantom divide cannot be realized. In particular, in the first three examples, the 
phantom crossing from the non-phantom phase to the phantom one and that from the 
phantom phase to the non-phantom one can occur. Moreover, in the last example, the 
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FIG. 21: wg as a function of t for / = 1, m = —2, n = —1, oq = —1/2, 6 = 1/2, d = 1 and /c = \/3. 



crossing of the phantom divide can be reahzed from the non-phantom (quintessence) phase 
to the phantom one (in which w < —1). Such a transition behavior is consistent with 
the observational suggestions 46145 1|. We also note that the qualitative behaviors of the 
numerical results shown in Figs. [T|-[2T] do not strongly depend on the initial conditions as 
well as the model parameters. Presumably, from such a nature it seems that even if we 
include matter, the matter-dominated stage would be realized and an attractor solution 
could be realized. 

We caution that the solutions reconstructed in this work might well be unstable. In 
addition, there appears divergence in the effective energy density in most of the examples 
during finite time. The stability of those models will be a crucial problem, which will spoil the 
models even if matter sectors are included. Accordingly, in order to solve this stability issue, 
we have to extend the present toy models to be more elaborate ones. Indeed, the resultant 
solutions could correspond to more complex models such as theories with non-linear kinetic 

n nn 

terms, e.g., the ghost condensate scenario [79|] and the Galileon gravity models |8CH84l|. It 
would be the most important future work of the present paper to examine the stability of the 
obtained solutions by comparing with the above theories which consist of non-linear kinetic 
terms, so that we can find some clue for the mechanism to make the derived cosmological 
solutions viable. 

Consequently, the divergence or oscillating behavior of effective energy density will dra- 
matically change the history of the universe, so that the process of the nucleosynthesis will 
probably be spoiled. Thus, it is crucial task for the present models to preserve the big 
bang nucleosynthesis (BBN), so that g-essence models can be realistic ones which are able 
to successfully resolve the mechanism of the current accelerated expansion of the universe. 
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Furthermore, we also note that the reconstructed expression of the potential V2 is given 
by a function of t and that of not (p. This is because it is difficult to obtain the explicit 
analytic form oi V2 = V2{ipip) and we find the analytic representation of and that of not 
(f). In this respect, the reconstructed procedure performed in this work would not perfect 
in order to derive the form of a g-essence model as a scalar field theory and investigate the 
effective mass through the potential expression. Since we concentrate on the cosmological 
evolution of the EoS for dark energy, the reconstruction of explicit form of the potential is 
an additional result of this work. In fact, it should be noted that what V2{ipip) is given as a 
solution of the equations of motion is very peculiar. In this work, we consider a toy model 
of g-essence in terms of the fermion ip phenomenologically. However, if we take a concrete 
fundamental theory in particle physics, it is considered that the particle physics theory can 
present the form of such kind of effective potential for ip. Furthermore, this situation would 
be similar to that in the quintessence models on the level that the potential cannot be given 
by particle physics because we do not have any particle physics theory which can give the 
very small mass of quintessence to realize the current cosmic acceleration, namely, in order 
for the current value of the quintessence potential to be equal to the very small value of the 
cosmological constant at the present universe. 

Finally, it is meaningful to emphasize the novel ingredients of this work. In g-essence 
models, there exist two dynamical components. One is the scalar field (p. Another is the 
fermion condensate u = ipi/j. These two components play a role of dark energy. In our 
analysis, the dynamics of both and u is included. As a consequence, the evolutions of 
cosmological quantities such as the deceleration parameter, the jerk parameter and the EoS 
in g-essence models become different form those in two scalar field models. This is one of 
the most important results obtained in this work. In addition, in Sec. IV B we explicitly 
illustrate that there exist g-essence models in which the crossing of the phantom divide can 
be realized. In general, from the quantum field theoretical view, in a single scalar field 
theory, if the EoS of the scalar field is less than —1, i.e., the universe is in the phantom 
phase, then the null energy condition is violated, so that the vacuum can be unstable, and 
hence there will appear microphysical ghosts. However, in g-essence models, in the field 
equation level there exist couplings of the fermion fields to H, and thus the dynamics of all 
the fields in the basic system is not exactly equivalent to that in a single scalar field theory. 
Accordingly, it is not so trivial whether the ghosts will exist in g-essence models examined 
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in Sec. IV B and the vacuum will be unstable. The investigations on the existence of the 
ghosts as well as the instabilities of the vacuum such as Laplace instabilities should be the 
future work in our studies of g-essence models with realizing the crossing of the phantom 
divide. The main purpose of this paper is to illustrate that the crossing of the phantom 
divide can occur in the framework of g-essence models, although there have been proposed 
other simpler models showing the crossing of the phantom divide in the literature. 

V. CONCLUSIONS 



In the present paper, we have studied the cosmological EoS for dark energy in g-essence 
models. In particular, we have found several g-essence models in which the universe is always 
in the non-phantom (quintessence) phase and thus the crossing of the phantom divide cannot 
occur. Furthermore, we have examined an explicit g-essence model in which the crossing 
of the phantom divide can be realized from the non-phantom phase to the phantom one. 
This transition manner is compatible with the analyses of the recent various cosmological 
observational data 
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m- 



We remark that in Ref . 6| , the limit on a constant EoS for dark energy has been analyzed 
as wde = —1.10 ± 0.14 (68% CL) in a flat universe. Moreover, in case of a time- dependent 

86| with wdeo 



EoS for dark energy obeying a linear form w-oE{a) = w^oEo + w^DEa (1 — a) 185 

and WDEa being the current value of wde and its derivative, respectively, the constraints 

have been estimated as wbeo = —0.93 ± 0.13 and wnK g. = —0-4:11l^'^1 (68% CL) |6|. If the 



more precise future CMB experiments such as PLANCK 87| support the phantom phase, the 
studies on the crossing of the phantom divide become more important in order to investigate 
the nature and property of dark energy. Thus, it can be considered that our results would 
be worthy of a clue to obtain the significant features of dark energy. 
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